Born Effective Charges and Infrared Response of LiBC 
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Calculations of the zone center optical mode frequencies (including LO-TO splitting) , Born effec- 
tive charges 7j* aa for each atom, dielectric constants eo and eoo, and the dielectric response in the 
infrared, using density functional linear response theory, are reported. Calculated Raman modes 
are in excellent agreement with experimental values (170 cm -1 and 1170 cm -1 ), while it will re- 
quire better experimental data to clarify the infrared active mode frequencies. The Born effective 
charges Z* a (i) have surprisingly different values for B and C, and (ii) show considerable anisotropy. 
Relationships between the effective charges and LO-TO splitting are discussed, and the predicted 
reflectivity in the range - 1400 cm -1 is presented. These results hold possible implications for Li 
removal in LiBC, and C substition for B in MgB2. 

PACS numbers: 77.22.-d, 78.30.-j, 63.20.-e 



I. INTRODUCTION 

Discovery of superconductivity near 40 K in 
and the rapid development in the under- 
standing of the microscopic mechanism (a specific 
type of electron-phonon coupling) has spurred ex- 
tensive work on its properties and the discovery of 
closely related materials. Interest in LiBC and its 
Li-deficient derivatives has been spurred by calcula- 
tional evidence @ that the electron-phonon coupling 
will be stronger, and consequently the superconduct- 
ing transition temperature is expected to be higher, 
than MgB2. LiBC is isostructural to, and isovalent 
with, MgB2 to the extent possible in a ternary com- 
pound, as we discuss more below. 

LiBC was first reported by Worle et al., 3] who 
provided the crystal structure, and also reported 
hole doping by removal of Li that changed it from 
reddish in color to black, and increased the con- 
ductivity. Recently there has been renewed inter- 
est in synthesizing this compound and studying its 
properties. Bharathi et al. reported synthesis for 
varying starting concentrations of Li,Q but did not 
report any samples becoming metallic. Hlinka and 
coworkers reported Raman scattering measurements 
on few-micron size crystallites, finding modes at 
170 cm -1 and 1176 cm -1 , which were recently con- 
firmed by Renker et al. (171 and 1167 cm -1 ).|(J 
Souptel et al. reported synthesis of samples for 
several flux concentrations, and like Bharathi et 
al. found conductivities characteristic of lightly 
doped semiconductors. Beginning optical studies 
on LiBC have been reported by Pronin et al. using 
an infrared microscope on small crystals, provid- 
ing infrared (IR) reflectivity in the 400-1900 cm -1 
range and measurements of the dielectric behavior at 
lower frequencies. We will return below to discussion 
of the implications of the Raman and IR data for the 
vibrational properties of LiBC. 



As for possible doping (de- intercalation of Li), all 
of these reports are probably consistent with very 
low doping, too low to drive the material conducting 
(in spite of the broad bands that promote doped- hole 
conduction). The susceptibility data indicate small 
concentrations of local moments, typical of semicon- 
ductors with point defects. Zhao, Klavins, and Liu, 
in their efforts at Li de-intercalation, observed evap- 
oration of Li upon vacuum annealing, with no sign 
of superconductivity in the annealed samples. 9] The 
amount of Li remaining in the sample has not yet 
been characterized, however (xrays are not sensitive 
enough) , so definitive results remain to be obtained. 



The purpose of this paper is to obtain in more 
detail the character of the intrinsic (insulating) 
phase of LiBC. Phonon frequencies calculated by 
linear response techniques using the full-potential 
linear muffin-tin orbital (FP-LMTO) method were 
reported by An et aL[lfJ In this paper we use re- 
lated by extended methods to calculate not only the 
phonon frequencies (again, with different codes) but 
the dynamic effective charge tensor Z* , the longitu- 
dinal - transverse (LO-TO) mode splittings, oscilla- 
tor strengths of the IR active phonons, and related 
dynamic dielectric behavior. Since Li is known, from 
calculation, to be ionic in this compound, one could 
anticipate Z* Li ~ +1, and because B and C lie side- 

Z* c and of 

-7i* Li from the acoustic sum rule. 



by-side in the periodic table, that Z B 



course Z* B + 7? c 



We find much more interestng behavior, however; 
Z* B and Z* c are of opposite sign, and large - they be- 
have like strongly charged cation and anion, respec- 
tively. We analyze this behavior, and also comment 
on the possible relevance this behavior has for the 
difficulty in de- intercalating Li from LiBC, and also 
for C doping in MgB2_ x C x . 




FIG. 1: Crystal structure of LiBC, showing the alternat- 
ing B-C arrangement in the graphene layers. B and C 
atoms also alternate along the c axis, resulting in two 
layers and two formula units per primitive cell. The 
Li ions reside in the interstitial positions between B-C 
hexagonal rings above and below. 



II. CALCULATIONAL METHODS 

The calculations were carried using the ABINIT 
code with Troullier-Martins pseudopotentials [111 
IT2I ]. Teter parametrization[l3j of the Ceperley- Alder 
exchange-correlation potential, and 216 k-points in 
the irreducible wedge of the Brillouin zone. The 
kinetic energy cutoff for the plane waves was 60 
Hartree. The /- sum rule (relative to unity) fl4| 
was satisfied with to 1.0005, which is one measure 
that the basis set is effectively complete. The ex- 
periment lattice constants (a = 2.752A, c = 7.058A) 
3] of LiBC (space group P63/ mmc, No. 194) were 
used. The calculation of effective charges follows the 
formalism of Gonze and Lee|l5| as implemented in 
the ABINIT code. 

The crystal structure of LiBC, pictured in Fig. 
is a direct generalization of that of MgB2. The B and 
C atoms form a flat graphene sheet, with B and C 
atoms alternating around each hexagonal unit. The 
stacking of B-C layers is alternating, such that each 
B has C neighbors along the c direction, and simi- 
larly for C. This stacking, which doubles the unit cell 
volume over that of MgB2, indicates that B-C bond- 
ing along c (to the extent that it occurs) is preferable 
to B-B and C-C bonding. Possibly this stacking re- 
sults in a favorable Madelung energy, although real 
static charges associated with B and C separately 
are almost impossible to define or calculate. The Li 
ions lie in the interstitial site (center of inversion) 
between the centers of the hexagonal units above 
and below, coordinated at equal distances with six 
B atoms and six C atoms. 



III. RESULTS 

A. Vibrational Frequencies 

The character of the various zone center modes 
in LiBC were presented earlier by An et [Toj 
For completeness, we repeat them here, with our 
calculated frequencies (in cm" 1 ), with macroscopic 
electric field contributions neglected; those correc- 
tions and longitudinal optic (LO) - transverse optic 
(TO) splittings are discussed below. Factor group 
analysis |l6j at the T point yields 15 optical modes : 
2^2« + 2i?ig + B^u with motions along the c direc- 
tion, and 2Ei u + 2E2g + E2 U with motions in the a — b 
plane (the latter type are all two-fold degenerate E 
modes). The 2A2 U and 2E\ U are infrared active and 
the 2E 2g modes are Raman active. [El Irl lei IToj 

• u=169: E2 S , B-C layers sliding against each 
other 

• w=299: B 2 g, B-C layers beating against each 
other along c 

• a>=292: E2U, Li layers sliding against each 
other 

• w=346: Ei u , Li layers sliding against the B-C 
layers 

• w=407: A2„, Li layers beating against the B-C 
layers, along c 

• w=510: Bi u , Li layers beating against each 
other along c 

• w=803: A2U, B-C puckering mode, all B atoms 
move oppositely to all C atoms, Li sites be- 
come inequivalent 

• w=829: B2 ff , B-C puckering mode, B moves 
with C atoms above/below it, Li sites remain 
equivalent 

• a>=1143: E itl , B-C bond stretching mode, the 
two layers are out-of-phase 

• o>=1153: E 2fl , B-C bond stretching mode, lay- 
ers in phase 

Our calculated values are also provided in Tabled 
where LO-TO splittings, where they arise, are also 
reported. Compared with the observed Raman ac- 
tive modes at 170 and 1176 cm _1 j^] and 171 and 
1167 cm -1 ,|(| the calculated values, 169 and 1153 
cm" 1 , are in excellent agreement. 

The LO-TO splitting for displacement parallel to 
the layers occurs for the "layer sliding" mode at 346 
cm" 1 (Auj = 25 cm" 1 ), and for the B-C "bond 



2 



TABLE I: Zone center optical modes of LiBC. The first 
column (E = 0) is for no electric field, the second is for 
the field lying in the plane (E _L c), and the third is for 
E || c. The 2Auu and 2E\ U are infrared active and the 
2E2 g modes are Raman active. Au> 2 — lo^q - cj^o 
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stretching" mode at 1143 cm -1 (Alo = 92 cm -1 ). 
For c polarization, the "layer beating" mode is at 407 
cm -1 (Aui — 100 cm -1 ), and the B-C "puckering" 
mode at 803 cm -1 is split by only Alu — 25 cm -1 . 
The large LO-TO splittings of the "layer beating" 
and "bond stretching" modes suggest they involve 
large effective charges, as we confirm below. 

Data on IR active modes have been published by 
Pronin et al. ||| who fit reflectivity data with modes 
at 540, 620, 700 (less clear) and 1180 cm -1 . The 
highest (and strong in the data) mode is close to our 
calculated value (TO 1143 cm -1 , LO 1235 cm -1 ), 
but the others are difficult to assign. The data, with 
three or four positions of structure, apparently in- 
clude contributions from c axis polarization, a — b 
plane polarization, and all angles between. Bharathi 
et al have presented IR absorption spectraQ with 
main peaks at ~380, 950, and 1200 cm -1 . The lower 
and upper of these are not far from our calculated 
values for E l9 modes (Table HJ, but the nearest IR 
mode to the 950 cm -1 peak is our A2„ mode at 803- 
828 cm" 1 . 



B. Effective Charges and Oscillator Strengths 

In a non-cubic lattice the Born (dynamical) effec- 
tive charge of an ion becomes a tensor. In hexago- 
nal symmetry such as in LiBC, the tensor is diago- 
nal and reduces to two values Z* x = Z* y = Z^ and 
Z* z = Z\. These charges, which show considerable 



anisotropy in LiBC, are given in Table II. The acous- 
tic sum rule J2 s ^s = is fulfilled to within 0.01, 
suggesting well converged calculations. In spite of 
the fact that they are neighbors in the periodic table 
and that both B and C are comfortable forming the 
graphene structure layer, when they alternate in this 
graphene layer their effective charges are vastly dif- 
ferent: substantial in size but different in sign. Given 
their similarities, the effective charge sum rule, and 
that Li is ionized, one might perhaps expect values 
close to +1, 0, and -1 (or possibly +1, and -^for 
Li, B, and C, respectively. They are, on (angular) 
average, much closer to the values +1, +2, and -3, 
which are not so far from the bizarre ionic configu- 
ration (only one with closed shells) of +1, +3, and 
-4 for Li, B, and C respectively that would reflect 
closed shells on all ions. It must be recognized, of 
course, that the Z* values are dynamic only, and re- 
flect the effects of covalency with respect to some 
reference ionic value; however, the reference values 
here are unclear (except for Li, which clearly should 
be +1). 

Due to the expected strong covalency in the B-C 
bonds, it perhaps is not surprising that the B and C 
effective charges are strongly altered from (smaller) 
reference ionic values. For c axis polarization, the 
Li value itself is unusual (almost 50% larger than its 
nominal value) , indicating Li is definitely involved in 
the interlayer coupling. We return to this important 
feature below. This partial Li covalency is proba- 
bly connected with the difficulty in deintercalating 
it from the LiBC lattice. [E 

Another view of the charge response can be ob- 
tained from the "mode effective charges" defined by 
Gonze and Lee^j|, which for the symmetry of LiBC 
can be written for each polarization as 

z* = y — „ . (i) 

Here each ion effective charge is weighted by the 
mode eigenvector that is normalized without the ion 
mass factor that occurs in the usual normalization 
of U Ka (k is the atom index, a is the Cartesian co- 
ordinate). The mode effective charge vectors indi- 
cate how strongly a mode will couple to an elec- 
tric field, and are listed in Table III for the IR ac- 
tive modes in LiBC. For the c axis polarized modes, 
|Z*| is 3.00 for the (softer) layer beating mode and 
2.10 for the (harder) B-C puckering mode. For pla- 
nar displacements, the (soft) layer sliding mode has 
the value 1.52, and the (hard) B-C bond stretching 
mode value is 5.50. These trends follow those found 
by Zhong et al. in perovskites.|l7| where displace- 
ments that modulated the "covalent" bonding pro- 
duced the largest mode effective charge. The value of 
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TABLE II: The Born effective charges of LiBC. The 
acoustic sum rule ^ s Z* = is fulfilled to within 0.01. 
= (1/3) TrZ* 





Z*(Li) Z*(B) Z*{C) 


zjj 
zi 


0.81 2.37 -3.17 
1.46 0.61 -2.07 


z* 


1.03 1.78 -2.80 



3.00 for the layer beating mode reflects the substan- 
tial Li "covalent" character for c axis displacements. 



C. Infrared Response 

The static electronic dielectric constants are cal- 
culated to be = 12.95, eL = 11.24. These values 
tend to be overestimated in LDA calculations such 
as those used here, presumably due to the under- 
estimation of the bandgap. Dielectric constants are 
also somewhat dependent (at the ~5% level) on the 
choice of pseudopotential. According to the gener- 
alized Lyddane-Sachs-Teller (LST) relation, 



n 



LO' ' 



J TO,m 



(2) 



which is used separately for each polarization, the 
static dielectric constants eo of both directions are 
17.4 and 18.5, respectively. In the GHz range (j~0.03 
cm -1 ), e s» 35 was reported by Pronin et aZ,|g but 
this value is not necessarily expected to be intrinsic. 

The relationship between the LO-TO splitting 
and the Born effective charges is given by the 
relational, (which holds for parallel and per- 
pendicular vibrations separately), 



47T 



J TO,m\ 



aa 



Vo 



E 



(Z* e) 2 
M K 



/ j ion,K ' 



(3) 



In this relation, m goes over the IR active modes of 
the given polarization, M K the ionic mass of aton re, 
Vo the volume of the primitive unit cell, and the right 
hand side has been expressed in terms of screened 
ionic plasma frequencies f^on.K that indicate how 
much each ion contributes to the combined strength 
of the resonances. For in-plane polarization, Li, B, 
and C contribute 6%, 26%, and 58%; the Li ion has 
minor effect and the response is dominated by B-C 
bond stretching. For c axis polarization, the con- 
tributions are 43%, 5% and 51% respectively; here 
the B is almost irrelevant, and the Li contribution is 
comparable to that of C. 



TABLE III: Mode effective charges of the IR active op- 
tical modes 



modes 


Directions 


a b c 


layer sliding Ei u 


0.65 -1.37 0.00 
1.37 0.65 0.00 


layer beating A2 U 


0.00 0.00 3.00 


B-C puckering Ai u 


0.00 0.00 -2.10 


bond stretching E^ u 


-5.50 -0.10 0.00 
-0.10 5.50 0.00 



The contributions can be examined mode by 
mode. The dielectric function can be expressed in 
terms of contributions from the IR active modes (m) 
as 



e(w) 



'TO.m 
,2 



= 1 



J TO,m 



J TO,m 



(4) 



The second line follows from the relation between 
the mode oscillator strength S m and the correspond- 
ing LO-TO splitting, provided in Table[I]in terms of 
Aw 2 . For in-layer polarization, large LO-TO split- 
ting of the bond stretching E\ u mode at 1143 (1235) 
cm -1 is the result of the large B and C effective 
charges; the splitting of the layer sliding mode at 
346 (371) cm -1 is an order of magnitude less. For c 
axis polarization, the effective charges are "derived" 
70% from the 407 (507) cm -1 layer beating mode 
with large LO-TO splitting involving the larger Z* 
for Li, and only 30% derives from the 803 (828) cm -1 
B-C puckering mode. 

From Eq. 0]the zero frequency limit is given by 



- 1 = 



- 1 



TO.m 



(5) 



which is different from the approximate "generalized 
Lyddane-Sachs-Teller" relation of Eq. [21 that is of- 
ten quoted. It is not as different as it looks. The 
"textbook" form of the dielectric function in terms 
of its zeroes and poles is 



n 



J LO,m 



W TO,m ~ ^ 



(6) 



While this expression has the same poles (at uJTO.m) 
as Eq. 21 the positions of its zeroes, at LOLO.m in 
this expression, are different from the zeroes in the 
general expression Eq. 01 However, zeroes have to 
occur between each of the poles, and the zeroes of 
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FIG. 2: Calculated reflectivity spectrum R(w) for LiBC 
for polarization in the a — b layer (top panel) and along 
the c axis (bottom panel), calculated from Eq. HU The 
dashed lines mark the positions of the TO and the LO 
modes for each IR active mode (two of them for each 
polarization). The damping was chosen to be 3% of the 
frequency for each mode. 

Eq. 0]are determined not only by uiLO.m but also by 
the polarization arising from other LO modes with 
ml 7^ to. The differences are small, however. For 
in-plane polarization, the zeroes of e(u>) from Eq. 
U] are 367 cm -1 and 1236 cm -1 compared to the 
calculate "LO" frequencies of 371 and 1235 cm -1 , 
respectively. (Note that the lower is pushed down, 
while the higher is increased, for the case of two LO 
modes as we have here.) For c axis polarization, Eq. 
0]gives 499 and 833 cm^ 1 compared to the calculated 
values of 507 and 828 cm -1 . These differences will 
be difficult to resolve experimentally. 



IV. DISCUSSION AND SUMMARY 



the spectra presented by Pronin et a/.|8| and by 
Bharathi et al.<$\ are not representative of pure po- 
larization data. 

Since it will be useful for comparison with single 
crystal data, in Fig. 2 we present our predicted re- 
flectivity R(uj) spectrum for both in-plane and per- 
pendicular polarization. The damping (which for a 
clean undoped sample would be due to anharmonic- 
ity) has been chosen for each mode to be 3% of the 
frequency (zero damping curves are not very repre- 
sentative). If the damping is not too large, and there 
is no reason to expect large anharmonicity in LiBC, 
then both TO and LO frequencies can be obtained 
with small uncertainty. 

It was mentioned in the Introduction that calcu- 
lations have indicated that, if partial Li removal can 
be achieved, Li^^BC should be a very good super- 
conductor. Such Li extraction is common in many 
materials, being the process that forms the basis for 
Li batteries. What our studies here have shown is 
that (primarily for c axis displacements) Li shows 
considerable covalency with the B-C layer, with Li- 
C coupling being the prominent feature. This par- 
tial Li covalency is probably connected with the dif- 
ficulty that several groups have foundf^, Q in de- 
intercalating it from the LiBC lattice. In addition, 
Cava, Zandbergen, and Inumaru reported briefly 
that LiBC seems to be highly resistant to chemical 
doping. 

The finding here that C is very different chem- 
ically from B in this system also carries some im- 
plications for C replacement of B in MgB2 . Recent 
studies confirm that about 10% of B can be replaced 
with C while retaining the structure, and the su- 
perconducting transition temperature remains high 
(T c =22 K).||a,|2l|22 If C alloying in the B sublat- 
tice were rigid band like, the additional 0.1 electron 
per cell would nearly fill the important a bands p3| 
that are responsible for superconductivity, and T c 
would be expected to vanish (or nearly so). The 
great difference between C and B behavior in LiBC 
reflects strong differences in bonding, and suggests 
something very different from rigid band "doping." 
Preliminary calculations for Mg B2-3iC x (ordered su- 
percells and virtual crystal) [24| indeed indicate sub- 
stantial non-rigid-band behavior: C is really differ- 
ent from "B with an extra electron." 



The calculated Raman-active frequencies are in 
excellent agreement with the data of two groups. 0, 
6] The IR modes remain to be verified and under- 
stood. The spectra are polarization dependent, and 



Although calculations such as those described here 
have been found in several systems to be accurate, it 
is highly desirable to obtain polarized single crystal 
IR data to confirm our predictions. 
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